Relaxation of antiferromagnetic order in spin- 1/2 chains following a quantum quench 
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We study the unitary time evolution of antiferromagnetic order in anisotropic Heisenberg chains that are 
initially prepared in a pure quantum state far from equilibrium. Our analysis indicates that the antiferromagnetic 
order imprinted in the initial state vanishes exponentially. Depending on the anisotropy parameter, oscillatory 
or non-oscillatory relaxation dynamics is observed. Furthermore, the corresponding relaxation time exhibits a 
minimum at the critical point, in contrast to the usual notion of critical slowing down, from which a maximum 
is expected. 



Introduction. Experiments with ultracold atoms offer a 
highly controlled environment for investigating open ques- 
tions of quantum magnetism. In particular, coherent spin dy- 
namics in a lattice of double wells has been observed in recent 
experiments, which have demonstrated remarkable precision 
in tuning magnetic exchange interactions [1]. The ability to 
observe quantum dynamics over long time intervals allows 
one to study strongly correlated states from a new perspective. 
The idea is to prepare the system in a simple quantum state 
which, in general, is not an eigenstate of the Hamiltonian, and 
investigate the dynamics that follows. In the two-spin sys- 
tem, studied in [1], the dynamics is completely tractable and 
describes simple oscillations between a singlet and a triplet 
states. 

In the present paper we investigate how the nature of the dy- 
namics changes in the case of a macroscopic number of spins 
interacting via nearest neighbor magnetic exchange. Are there 
new effects, and in particular new time scales, dynamically 
generated by the complex many-body evolution? Our starting 
point for investigating this question is the spin-^ anisotropic 
Heisenberg (or XXZ) model on a one-dimensional lattice 

j 

This model provides a good effective description of two- 
component Bose or Fermi systems deep in the Mott insu- 
lating phase. The interaction parameters are dynamically 
tunable [2], realizing ferro- (J < 0) or antiferromagnetic 
(J > 0) couplings over large ranges of the anisotropy param- 
eter A > 0. We take the initial state to be a perfect anti- 
ferromagnetic (Neel) state iV'o) = I Tit ••• iti)- Such a 
state has been achieved with high fidelity by Trotzky et al. 
[1] using decoupled double wells. Note that | V'o ) is the 
ground state of the Hamiltonian with A = oo. We study the 
subsequent time evolution of the staggered magnetic moment 
ms{t) = ^Ej(-lP(V'o|5'|W|V'o) under the influence of 
the Hamiltonian (1) at different values of anisotropy A using 
a numerical matrix-product method [3]. The dynamics is in- 
dependent of the sign of J and the results are valid for both 
ferro- and antiferromagnetic couplings. To substantiate our 
findings, we consider another, closely related model, given by 



the XZ-Hamiltonian [see Eq. (5)], which allows exact calcu- 
lation of the dynamics and displays similar behavior. 

Theoretical interest in this class of problems, known as 
quantum quenches [4, 5, 6, 7, 8, 9], has been invigorated by 
advances in experiments with ultracold atoms [10]. In par- 
ticular, macroscopic order parameter oscillations have been 
predicted to occur following a quantum quench in a variety 
of such systems [11, 12, 13, 14]. We shall see that such os- 
cillations are also found in the XXZ-chain with easy-plane 
anisotropy (A < 1), and that they are essentially the same 
as the singlet-triplet oscillations observed in the two-spin sys- 
tem [1]. Accordingly, the oscillation frequency is directly re- 
lated to the magnetic exchange interaction J. More impor- 
tantly, for non-zero A we find a fundamentally new mode of 
many-body dynamics which always leads to exponential de- 
cay of the staggered moment regardless of whether the short- 
time dynamics is oscillatory or not. In contrast with the oscil- 
lation frequency, the relaxation time is an emergent scale gen- 
erated by the highly correlated dynamics and hence cannot be 
simply related to the microscopic parameters. We find a di- 
verging relaxation time in the two limits A ^ and A ^ oo. 
Of particular interest is the relaxation time at the isotropic 
point A = 1, which for the ground state properties marks 
a quantum phase transition from a gapless "Luttinger liquid" 
phase (A < 1) to a gapped, Ising-ordered antiferromagnetic 
phase (A > 1). Interestingly, the relaxation time is minimal 
in the vicinity of the critical point, where its value is simply 
determined by the magnetic exchange interaction r ~ 1/ J. 
This accelerated relaxation stands in remarkable contrast to 
the notion of critical slowing down, valid for a small pertur- 
bation of the order parameter from equilibrium. In fact, if the 
prepared initial state is close to the equilibrium state, then the 
relaxation time of the order parameter is expected to diverge 
as the system approaches the critical point [15]. We find an 
opposite trend in the dynamics of the prepared Neel state. In 
the long-time limit, our results suggest that local magnetic or- 
der vanishes for all values A < oo. 

The solution of the quench dynamics in the XXZ-chain 
involves in principle all energy scales of the Hamiltonian 
and approximative methods become essentially inaccurate in 
many cases. The mean field approximation for example leads 



2 4 6 



10 12 14 16 18 




FIG. 1 : Dynamics of the staggered magnetization {t)m the XXZ- 
chiain (a) and the XZ-chain (b). Symbols correspond to numerical 
results, lines represent analytical results or fits by corresponding laws 
(see text). For A = the typical behavior of the error is illustrated 
by comparing the numerical iTEBD result with 2400 retained states 
to the exact curve: the absolute deviation from the exact curve is less 
than 10"" for t < tmnaway For A 7^ data beyond tmnaway is 
omitted. 



to contradictions with our results - an algebraic decay for 
A < 1 and a non-vanishing asymptotic value of the staggered 
moment for A > 1 [14]. Renormalization group based ap- 
proaches [5] are restricted to low-lying modes, which is not 
sufficient in the present case. The exact numerical results pre- 
sented in this study go further than the predictions of low- 
energy theories [5]. 

Before delving into a more detailed study it is instructive 
to consider the so-called XX-limit (A = 0) of the Heisen- 
berg chain ( 1 ), which can be mapped onto the problem of free 
fermions. In this case one easily obtains an analytic expres- 
sion for the time evolution of the staggered magnetization: 
ms{t) = Jo{2Jt)/2 (Fig. 1). Here Jo denotes the zeroth- 
order Bessel function of the first kind. Thus, after a short 
transient time t ^ J^^, the staggered magnetization displays 
algebraically decaying oscillations originating from the finite 
bandwidth of the free-fermionic model, 



^ : cos(2 Ji - -) . 

/4TTt 4^ 



(2) 



In general, we are interested in generic behavior of the relax- 
ation dynamics on large time scales. We adopt a definition of 
relaxation which does not rely on time-averaged equilibration 
of the observable, but instead requires exact convergence to 
the asymptotic value, as defined in Ref. [16]. From this point 
of view, the oscillations in the XX-Iimit are characterized by 
an infinite relaxation time. 

XXZ-model. In the general case of A 7^ 0, the prob- 
lem is no longer analytically treatable and we have to re- 
sort to numerical techniques. We use the infinite-size time- 
evolving block decimation (iTEBD) algorithm [3], which 




FIG. 2: Absolute value of the staggered magnetization in the XXZ- 
model. Symbols represent numerical results, solid curves correspond 
to fits by the exponential law (3), straight lines point out the expo- 
nential decay. 
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FIG. 3: (a) Focus on the XXZ-chain close to the critical point A — 
1. (b) Comparison of the XXZ-chain (symbols) and the XZ-chain 
(dashed lines) for strong anisotropics, solid lines correspond to an 
exponential fit. The dynamics of the staggered magnetization of the 
XXZ- and XZ-chains converge towards each other in the large-A 
limit. 



implements the ideas of the density matrix renormalization 
group (DMRG) method [17] for an infinite system. The al- 
gorithm uses an optimal matrix-product representation of the 
infinitely extended chain, keeping only the dominant eigen- 
states of the density matrix of a semi-infinite subsystem, in 
combination with a Suzuki-Trotter decomposition of the evo- 
lution operator. This method is very efficient for small t, how- 
ever, the increasing entanglement under time evolution [18] 
requires to retain an exponentially growing number of eigen- 
states. We find that the error of our calculations behaves in a 
similar way to that of the finite-size DMRG algorithm and the 
methodology developed in Ref. [19] can be applied in order 
to control the accuracy [20]. By carefully estimating the run- 
away time via comparing results with different control param- 
eters [19], the absolute error in the plotted data is kept below 
10"®. Using 2000 states and a second-order Suzuki-Trotter 
decomposition with a time step 5 ~ 10"'^ for large A and 
up to 7000 states with 5 ~ 10"^ J"^ for small A, an interme- 
diate time regime Ji < 16 can be reached, which in general 
far exceeds the short transient time. 

An overview of the results is presented in Fig. la. For small 
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anisotropics we find oscillations of the order parameter similar 
to those in the XX-limit, but with a decay time decreasing 
upon approaching the isotropic point A = 1. In the easy- 
axis regime A > 1 of the XXZ-model, the relaxation slows 
down again for increasing A, and we observe non-oscillatory 
behavior for A ^ 1. 

Fig. 2 focuses on easy-plane anisotropy < A < 1. The 
results for < A < 0.4 are well described, for accessible 
time scales, by exponentially decaying oscillations 

nisit) (X e'*^'' cos{ujt + <j>) . (3) 

The oscillation frequency is almost independent of the 
anisotropy, while the relaxation time r increases with decreas- 
ing A. Logarithmic divergence of the relaxation time in the 
limit A ^ is suggested by the fit shown in Fig. 4a. The 
picture is less clear closer to the isotropic point. For the range 
0.5 < A < 1 there appears to be an additional time scale 
after which the oscillations start to decay even faster than ex- 
ponentially, simultaneously the period of the oscillations is 
reduced. Therefore, the relaxation times plotted in Fig. 4a are 
only valid within an intermediate time window, whose width 
shrinks upon approaching the critical point. 

For intermediate easy-axis anisotropics 1 < A < 3, the 
magnetization does not reach a stable regime within the nu- 
merically accessible time window (Fig. 3a). The compUcated 
behavior of ras{t) in this parameter range can be ascribed 
to the interplay of processes at all energy scales. Neverthe- 
less, the numerical data suggest that the relaxation is fastest 
close to the isotropic point, in the range between A = 1 and 
A = 1.6. A simple generic type of behavior is recovered for 
large anisotropics A > 3. The numerical data in Fig. 3b indi- 
cates exponential relaxation of the staggered magnetization 

m,(t) cx e~*/^. (4) 

The relaxation time scales roughly quadratically with A (Fig. 
4a). Oscillations do persist on top of the exponential decay, 
but they fade out quickly. 

XZ-model. We now turn to the study of the XZ- 
Hamiltonian, 

i^xz = J ^ { 2 + ^ ^l^l+i } • (5) 

3 

In this model a quantum phase transition separates two gapped 
phases at Ac = 2, with antiferromagnetic order in z-direction 
for A > Ac and in x-direction for A < Ac. Unlike the 
XXZ-model, it can be easily diagonalized analytically. In or- 
der to study the staggered magnetization of the XZ-model, we 
have to calculate the two-spin correlation function C(n, t) = 
(— l)"('0o|5'o (i)5'^(t)|'0o) in the infinite-range limit, since 
mg{t) = liuin^oo C{n,t). Using standard techniques (see 
[4] and references therein), we express this two-spin correla- 
tor as a Pfaffian, with coefficients calculated in a similar man- 
ner as for the Ising model in a transverse field [13]. Exploiting 
the light-cone effect [18, 21], we are able to evaluate numer- 
ically the order parameter dynamics up to times of the order 
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FIG. 4: Relaxation time r and oscillation period T — ^ as a func- 
tion of anisotropy in the XXZ- and XZ-models. Logarithmic or al- 
gebraic laws are emphasized by solid lines. In the region close to the 
critical point of the XXZ-model (indicated by the question mark) it 
becomes impossible to extract a relaxation time from the numerical 
results. 

of Jt Ri 100. The results are displayed in Fig. lb. An ana- 
lytic expression can be derived for A = 0, which is given by 
ms{t) = 0.5cos^(Ji). For A < Ac, exponentially decaying 
oscillations 

ms{t) oc e^*/'^(cos^(ijji) — const.) (6) 

reproduce the numerical results at large times very well. For 
A > Ac, the staggered magnetization decays exponentially 
with no oscillations at large times [Eq. (4)]. In contrast to the 
XXZ-model, the oscillation period in the XZ-model diverges 
at the isotropic point A = Ac and the latter exactly marks the 
crossover between oscillatory and non-oscillatory behavior of 
TOs(t). We have extracted the relaxation times from exponen- 
tial fits to the numerical data, showing a clearly pronounced 
minimum right at the isotropic point (see Fig. 4b). The relax- 
ation time scales as r cx A^^ for A < Ac and as r cx A^ for 
A > Ac. 

Apart from the numerical evaluation of the Pfaffian, we 
can prove rigorously that in the infinite-time limit the stag- 
gered magnetization vanishes for all anisotropics in the range 
Ac < A < oo. Indeed, since the Pfaffian reduces to aToeplitz 
determinant at t ^ oo [13], we can use Szego's lemma to 
calculate the large-distance asymptotics of the two-spin cor- 
relator in the above-mentioned regime, obtaining for n ^ 1, 

limt^oo C{n, t) ^ J ^ i-i-\/i-4 /A_^ ^ which immediately 

implies that ms{t — > oo) — 0. 

Discussion. We have analyzed the dynamics of the stag- 
gered magnetization in the XXZ- and XZ-models following a 
quantum quench. Our main result is that in both models there 
is a dynamically generated relaxation rate which is fastest 
close to the critical point. This point also marks a crossover 
between oscillatory and non-oscillatory dynamics of the or- 
der parameter. The dynamics of the magnetic order parameter 
turns out to be a good observable for the quantitative extrac- 
tion of non-trivial time scales. In general, this is not possible 
from other observables such as correlation functions, which 
reveal interesting features such as the horizon effect [5] but ex- 
hibit only slow relaxation dynamics [8, 20]. Furthermore, we 
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have focused on the Neel state as an experimentally relevant 
initial condition. We point out however, that our results are 
generic and hold for all antiferromagnetic initial states with 
sufficiently small correlation length [20]. 

The existence of a minimal relaxation time at the critical 
point is opposite to what one would expect from the phe- 
nomenon of critical slowing down of order parameter dy- 
namics near equilibrium. In the XZ-model and the easy-axis 
phase of the XXZ-model, where the excitation spectrum is 
gapped, the effect can be understood using a phase-space ar- 
gument: the relaxation of the initial state is dominated by scat- 
tering events between high-energy excitations introduced into 
the system through the initial state. As a result of the exis- 
tence of the gap, the phase space for scattering events is re- 
stricted. This leads to an increasing relaxation time as the gap 
increases, whereas one expects a minimal relaxation time at 
the critical point, where the gap vanishes. The above argu- 
ment can not be applied directly to a quench into the gapless 
phase in the easy-plane regime of the XXZ-model. Rather, the 
situation seems to be similar to a quantum quench of the Bose- 
Hubbard model from a Mott insulator to a superfluid phase. In 
the latter case, oscillations of the superfluid order parameter 
have been predicted, with a damping rate that diverges at the 
critical point in one and two dimensions [11]. 

The absence of a sharp signature of the quantum phase tran- 
sition in the XXZ-chain prepared in a Neel state is in contrast 
with what one has in the case of the initially prepared ferro- 
magnet with a single kink-impurity, studied for example in 
Ref. [19], where the two phases are characterized by clearly 
distinct transport properties. We note that this initial state is 
much closer to the ground state of the Hamiltonian and the 
important energy scales are considerably smaller than in the 
case of the initial Neel state. The opening of an exponentially 
small gap at the phase transition is therefore more likely to be 
relevant. 

The time evolution of an initial state which is equivalent to 
the Neel state has been recently studied by Cramer et al. [7] in 
the context of the one-dimensional Bose-Hubbard (BH) model 
with on-site repulsion U as interaction parameter (the equiva- 
lence becomes apparent in the fermionic representation of the 
XXZ-Hamiltonian). Although the BH Hamiltonian itself and 
the XXZ-model share some properties in the non-interacting 
limit, there is one substantial difference: in the BH model at 
half-filUng no equilibrium critical point is crossed by chang- 
ing the interaction U and the symmetry-broken initial state 
never becomes the ground state. Unlike the XXZ-model, the 
oscillations of the local observable in the BH model appear to 
be decaying algebraically for all values of interaction U and 
no crossover to a non-oscillatory regime has been observed 
[7]. These differences point out the crucial role of the equi- 
librium phase transition to the reported behavior of the order 
parameter dynamics of the XXZ-chain. 



Experimental results [1, 22] suggest that effects of density 
fluctuations beyond second-order magnetic exchange may be 
important for reproducing the dynamics in full detail. This 
statement is also supported by very recent numerical results 
[23]. Nevertheless we expect that our main result, the exis- 
tence of a minimum in the dynamically generated relaxation 
time close to the critical point, is insensitive to these details. 
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